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Abstract. We investigate the influence of different magnetic field configurations on the temperature distribution 
in neutron star crusts. We consider axisymmetric dipolar fields which are either restricted to the stellar crust, 
"crustal fields" , or allowed to penetrate the core, "core fields" . By integrating the two-dimensional heat transport 
equation in the crust, taking into account the classical (Larmor) anisotropy of the heat conductivity, we obtain 
the crustal temperature distribution, assuming an isothermal core. Including quantum magnetic field effects in 
the envelope as a boundary condition, we deduce the corresponding surface temperature distributions. We find 
that core fields result in practically isothermal crusts unless the surface field strength is well above 10^^ G while 
for crustal fields with surface strength above a few times 10^^ G significant deviations from isothermality occur 
at core temperatures inferior or equal to 10* K. At the stellar surface, the cold equatorial region produced by the 
quantum suppression of heat transport perpendicular to the field in the envelope, present for both core and crustal 
fields, is significantly extended by the classical suppression at higher densities in the case of crustal fields. This can 
result, for crustal fields, in two small warm polar regions which will have observational consequences: the neutron 
star has a small effective thermally emitting area and the X-ray pulse profiles are expected to have a distinctively 
different shape compared to the case of a neutron star with a core field. These features, when compared with 
X-ray data on thermal emission of young cooling neutron stars, will open a way to provide observational evidence 
in favor, or against, the two radically different configurations of crustal or core magnetic fields. 
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■""j 1. Introduction regions of the core where neutrons are superfluid and pro- 

r> I ton superconducting. 
■ The presence of strong magnetic fields in neutron stars 

" " ' is one of their distinctive characteristics. In typical neu- Any observed magnetic field, either based on crustal or 
tron stars the observed and/or inferred surface fields are core currents, has to penetrate the crust matter, thereby 
of the order of 10^^- -^^ G ; for magnetars they even reach affecting its transport properties. Moreover, the struc- 
2^014. ..15 Q rp^^^ strength, structure and time evolution of ture of the field in the shallow spherical shell layer below 
the magnetic field is intimately related to its origin, which the surface has a direct effect on the surface tempera- 
is still an open problem. Basically, two qualitatively dif- ture and leads to a non-uniform distribution with obser- 
ferent types of field structures are conceivable, one having vational consequences (Page 1995) if its strength is above 
an initial field penetrating the whole star while the other 10^" G. Roughly, the effects of the magnetic field onto 
is characterized by having the field and its supporting cur- the transport processes can be divided into classical and 
rents restricted to the stellar crust (see, e.g., Chanmugam quantum ones (see, e.g., Yakovlev & Kaminker 1994 for 
1992). To date, there is still no compelling observational a review). The classical effects are due to the Larmor ro- 
evidence in favor of or against either of these two hypothe- tation of the electrons, the main carriers of charge and 
ses but recently Link (2003) argued that long period pul- heat, and are determined by the magnetization parame- 
sar precession, as observed in PSR B1828-11 (Stairs et al. ter uJbT where uJb = eB/mlc is the gyrofrequency of the 
2000), may be impossible if the magnetic field penetrates electrons, r being their relaxation time and m* their ef- 
fective mass. Quantized motion of the electrons transverse 
to the magnetic field causes the quantum effects which are 



Send offprint requests to: U. Geppert, e-mail: urme@aip.de of importance only if few Landau levels are occupied, re- 



2 



Geppert, Kiiker & Page: Temperature distribution in magnetized neutron star crusts 



quiring thus densities p < pB = 7045(5/10^^0)3/2 g 
and temperatures T <^ Tb ~ 1.34 • 10^(B/10^2g) 
K (Potckhin ct al. 2003) where < A> and < Z > are the 
average mass and charge numbers of the ions. From these 
numbers it becomes clear that quantum effects will play 
an important role for strong magnetic fields in the outer- 
most layers of the neutron star crust - the thin low density 
shell of the envelope - but are negligible in the deeper lay- 
ers. The thinness of this envelope justifies a study of heat 
transport in a plane parallel, one dimensional approxima- 
tion and many such calculations have been performed (see, 
e.g., for the most recent one, Potekhin & Yakovlev 2001, 
hereafter PYOl, and references therein). Two-dimensional 
calculations of heat transport with magnetic field have 
been presented by Schaaf (1990a,1990b) who however re- 
stricted himself to the thin envelope and an uniform mag- 
netic field. Tsuruta (1998) has presented results of two- 
dimensional cooling calculations of neutron stars which 
included the quantizing effect of a dipolar magnetic field 
in the envelope. These 2D calculations showed that the ID 
approximation is indeed very good when the field aflfects 
heat transport only in the thin envelope. 



The deeper layers of the crust can, however, be af- 
fected by classical effects in case of strong fields, i.e., large 
uJb, and/or low temperature, i.e., high r. In such a case 
the usual assumption of an isothermal crust is question- 
able and this is the issue we want to address in this pa- 
per. In particular, in the case of a crustal magnetic field 
its strength in the crust unavoidably exceeds its surface 
value by one to two orders of magnitude (see, e.g.. Page 
et al. 2000) and its effect can naturally be expected to 
be stronger than in the case of a field permeating the 
whole star. This may open a way to study the structure of 
the magnetic field in the crust and provide observational 
features which will allow us to discriminate between the 
above mentioned two types of hypothesized field structure, 
crustal vs. core. 



The paper is organized as follows: In the next Sect. 
2 the basic equations are introduced which describe the 
magnetic field and the heat transport infiuenced by it. The 
components of the heat conductivity tensor are given and 
the outer boundary condition is discussed. The physical in- 
put as well as the numerical method are shortly described. 
In Sect. 3 the results of the numerical calculations are 
presented. For different core temperatures, magnetic field 
strengths and geometries the crustal temperature profiles, 
the surface temperature distributions and the correspond- 
ing luminosities are calculated. Sect. 4 is devoted to the 
discussion of the consequences of the magnetic field effects 
on to the crustal and surface temperature distribution. 



2. Physics Input 

2.1. Heat transport and magnetic field 

The thermal evolution of the crust is determined by the 
energy balance equation 



dT 

C — = Sources - Sinks - V • F 
ot 



and the heat transport equation 



-k-VT 



(1) 



(2) 



where C denotes the specific heat per unit volume, F is 
the heat fiux density and k the tensor of heat conductivity. 

In this paper we intend to consider only the effect of 
the crustal magnetic field on to the stationary temperature 
distribution in the crust, which, in a first approximation, is 
assumed to be free of heat sources and sinks. The cooling 
process itself as well as the back reaction of the now non- 
spherically symmetric temperature distribution on to the 
magnetic field decay will be beyond the scope of this work. 

In the relaxation time approximation the components 
of k parallel and perpendicular to the magnetic field, K|| 
and K± respectively, as well as the Hall component, ka, 
are related to the scalar heat conductivity kq and to the 
magnetization parameter cJbT by (Yakovlev & Kaminker 
1994) 



For Kq we use 



1 + {uJnT)^ 



Kq = 



Ka = LObTK± (3) 



(4) 



where u = 1/r is the effective electron collisional fre 
quency and is given by the sum 



(5) 



where i/ph is the effective collisional frequency for electron- 
phonon, fion for electron-ion, and Vi^p for electron- 
impurity collisions. Since for the range of densities and 
temperatures we are interested in the crust is always in 
the solid state we neglect Uion and for z/ph we use the calcu- 
lations Baiko & Yakovlev (1996) while we take t'imp from 
Yakovlev & Urpin (1980). Everywhere in this work we 
assume an "impurity parameter" Q = Xi^p{AZy equal 
to 0.1, where Ximp is the fractional number of impuri- 
ties which have a mean square charge excess of {AZy. 
Figure 1 shows the resulting kq: note that Vimp is almost 
temperature independent while Vph scale approximately as 

so that if = fph then the "phonon-only conductivity" 
1^0 ph ^ (Fig. 1 left panel) but if = i^^^p we would 
have an "impurity-only conductivity" Kq imp ^ (Fig. 1 
right panel) and the exact Kq (Fig. 1 central panel), which 
one could write as kq = {I/kq ph + l/Koimp)~^, shows both 
types of behaviours depending on which type of collision 
process dominates. 

The magnetization parameter iVbT varies strongly 
throughout the crust by many orders of magnitude. In 
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Fig. 1. Thermal conductivity /to, in cgs units, versus density p and temperature T. The left panel shows the phonon-only 
contribution /to ph, the right one the impurity-only contribution /to imp, and the central panel the complete /to (with an impurity 
concentration Qimp = 0.1). 



lOb both B and m* span a large range of values through- 
out the crust and in r there is also a strong dependence on 
the temperature T, chemical composition and the thermo- 
dynamic phase of the matter. For illustration we show, in 
Fig. 2, WbT in the crust at various uniform temperatures 
and a uniform, field strength: however, neither T nor B 
will be uniform in our realistic calculation presented be- 
low. Note that values of at T = 10^ and 10^ K are 
very close to each other because at such low T is dom- 
inated by fc'iinp which is temperature independent while 
when going to increasingly higher temperatures i^ph con- 
tributes more and more and hence r decreases. Finally, 
Fig. 3 shows the resulting /tj_': notice its very different T- 
p behaviour compared to kq due to the simple fact that 
at high WbT we obtain ki_ oc while /to oc r. 

2.2. The Crustal Magnetic Field 

A dipolar poloidal magnetic field can be conveniently de- 
scribed in terms of the (possibly time dependent) Stoke's 
stream function S = S{r,t). The vector potential A = 
(0, 0, A^) is written as A^ = S{r, t) sm9/r, where r, 9 and 
(f are spherical coordinates. The field components are then 
expressed as 



Br = + 



2 008 6* 



Sir,t) = Bo 



cos 9 



s{x,t) 



Ba 



sine dS{r, t) 



dr 



Bq sva.6 ds{x,t) 
2 X dx 



B^ = 



(6) 

(7) 

(8) 

where Bq is the field strength at the magnetic pole, and 
we have introduced normalized variables 5 = sBoR%g/2, 
with ,s = 1 at the magnetic north pole, and x = r/R^qs- 
The vacuum solution, outside the star, is simply s = 1/x. 
At the surface, x = 1, the standard boundary condition 
for matching to a vacuum dipole field should be fulfilled, 

ds_ „ . ds 

dr Rns 







I.e. 



dx 



(9) 



and at the center regularity requires S{r = 0,t) =0. 

For the present investigation we select for the crustal 
field configuration a snapshot of the evolution of S{x,t). 
The Stoke's function S{x, t) was calculated by solving the 
induction equation, applying the above boundary condi- 
tions and an electric conductivity which reflects the same 
microphysics as the heat conductivity does for the model 
under consideration (see, e.g., Geppert fc Urpin 1994 or 
Page et al. 2000). 

The initial value S{x, t = 0) is a priori unknown and 
depends on the field generation process. For a crustal mag- 
netic field, the function S{x, t = 0) initially vanishes in the 
core and, due to proton superconductivity, the Meissner- 
Ochsenfeld effect prevents the field from penetrating into 
the core (sec, e.g. Page et al. 2000). A typical crustal field 
structure is shown in the right panel of Fig. 4. In case the 
magnetic field is maintained by electric currents circulat- 
ing in the core the field penetrates the crust too but has 
a qualitatively different structure. Let us assume that for 
the core field there are no currents in the crust and the 
field is maintained in the core by axisymmetric currents 
circulating around the center of the star. Then a dipolar 
field will penetrate the crust with the components 



= Bc 



cos 6 



Bq sin 9 




(10) 



(11) 

(12) 



where Bq denotes again the polar surface value of the 
magnetic field. The crustal field lines of such a field are 
prc;scnted in the left panel of Fig. 4. A comparison of 
the geometries of both fields immediately suggests that 
the star-centered core field will not significantly affect the 
heat transport through the crust while a crustal field may 
cause drastic changes. Moreover, in the crust Be has to 
be much larger than B'g for a given identital external 
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field, since the flux of a crustal field is compressed into a 
layer of thicknes Ar — ^ of iJjvs, while the flux of a 
core field can expand within the whole core, i.e., typically 



2.3. The Two-Dimensional Heat Transport Equation 

Denoting the components of the temperature gradient VT 

parallel and perpendicular to the unit vector of the mag- 
netic field, b as well as the Hall, component by 

(VT)|| =6(VT-6) , (VT)i = 6 X (VT X 6) , (13) 
(Vr)A = b X VT , 

the magnetic-field-influenced heat flux can be expressed 

by 

k-VT = 

K\\ b (VT ■b) + K±bx (VT x 6) + ka & x VT . (14) 

This provides the following expression for the heat flux in 

terms of the scalar heat conductivity, the magnetization 
parameter, the temperature gradient and the unit vector 
of the magnetic field: 



its radial derivatives (see Eqs. 6, 7) and by introducing 
the heat conductivity coefficients 



Xi 



1 + (a;«T)2 ' ^2 1 + (u;,r)2 ' 

the radial and meridional components of the magnetic 
fleld-dependent heat flux have the following form: 

dT 



(16) 



dT s ds 



dT s2 



ox 



(17) 



idT 

X oO 



dT 



ds 



dT 1 



d6 4a;3 V dx 



ds 



sin^ e , (18) 



F = -K • VT : 



1 + {oJsr? 
[VT + [ustY b (VT ■ 6) + 6 X VT] 



(15) 



The last term on the r.h.s. represents the Hall component 
of the heat flux. Its divergence vanishes as long as the mag- 
netic field and the temperature gradient are assumed to 
be axially symmetric, i.e., do not depend on the azimuthal 
angle tp. By use of the representation of the dipolar mag- 
netic field in terms of the Stoke's stream function s and 

1000 F 



where T is measured in units of Tcorc and the heat flux is 
normalized on T^ore/ Rns- Note that in Eq. 16 the magne- 
tization parameter uJe^t is calculated by use of the polar 
surface field strength B^. 

Axial symmetry implies that, along the polar axis, Fg 
and dFg/d9 vanish while x^F^ is constant. Mirror sym- 
metry at the magnetic equator implies that Fg vanishes 
(but dFo/dO can be very large). 

Our aim in this paper will be to find stationary solu- 
tions for the temperature distribution in the neutron star 
crust, i.e. to solve the equation 



V F 



1 d{x^F^) 1 d{sm.eFe) 



dx 



+ 



XSIQ.6 



d6 



= 



(19) 




Density [g/cm"'] 



Fig. 2. Magnetization parameter ujbt vs. density at six diflter- 
ent temperatures (as labeled on the curves) assuming a uniform 
magnetic field of strength B = 10^^ G. Its value for different 
field strengths scales linearly in B. 



7 LogT(K) 



Fig. 3. Thermal conductivity, in cgs units, perpendicular to the 
magnetic field, kx, vs. density p and temperature T for a uni- 
form magnetic field of strength 3 x 10^^ G. For field strengths 



» 10 G, i.e., for lobt » 1, kx scales as B 
parameter Q is assumed to be 0.1. 



The impurity 
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Fig. 4. Magnetic field lines of the core field (left panel) and of 
the crustal field as applied in our calculations. Both field con- 
figurations match for r > Rns with a dipolar poloidal vacuum 
field. They are presented in a meridional plane, of the crust for 
Pn > p > Pb (see Sects. 2.4 and 2.5). The polar surface value 
Bo is identical for both field types. 

The boundary condition for this equation at the crust-core 
boundary is " T fixed and uniform " and at the surface it 
as discussed in the following subsection. 

2.4. Outer boundary: Ti,{B) - Ts{B) relationship 

In the lowest-density layers, close to the surface, matter 
is no longer degenerate and the magnetic field affects the 
equation of state. Appropriate treatment of this layer re- 
quires solving for hydrostatic equilibrium simultaneously 
with heat transport. To avoid this problem, we separate 
this layer, cafied envelope, from the crust and incorporate 
it in the outer boimdary condition. We stop the integra- 
tion at an outer boundary density ph, at radius rt,, and, for 
each latitude 9, obtain a radial flux Fr{pb,0) and a tem- 
perature T(pi,, 9) which we match to the envelope. In the 
envelope approximation, a surface temperature Tg, and 
hence an out coming fiux F"^^ = crTf, is chosen and hy- 
drostatic equilibrium and heat transport are solved toward 
increasing densities up to pb giving the temperature T™^ 
at that density. Varying Tg gives a "Tf, — Tg relationship" . 
The matching of the envelope with our interior calcula- 
tion is simply obtained by imposing T™^ = T{ph,9) and 
^env ^ ^rp4: _ p^(^p^^ which is our outer boimdary con- 
dition. The magnetic field strongly affects the structure 
and transport properties of the envelope and the "T^ — Tg 
relationship" : at each point {rt , 9) we apply an envelope 
with a magnetic field B equal to the field we have at that 
point, hence giving us a Tb{B) — Ts{B) relationship. 

Many calculations of magnetized envelopes have been 
presented and we use the latest and most accurate one 
from PYOl. These authors pose the bottom of the enve- 



lope at the neutron drip density, i.e., ~ 4 • 10^^ g cm~^, 
while we intend to apply lower values down to pb = 10^° 
g cm~^ in order to extend the 2-dimensional transport 
calculation as far as possible while still being safely in 
the non-quantizing regime. Since the temperature profile 
in the 1-dimensional envelope calculations of Potekhin & 
Yakovlev is quite flat between 4 • 10^^ g cm~^ and 10^" 
g cm~^ the same Tb — Tg relationship can be applied for 
the lower pb which we consider a good approximation. 
Explicitly, we apply their equations (26) to (30) but re- 
place their Tcore which they assume to be spherically sym- 
metric by our calculated angle-dependent Tb{9). 

To illustrate the main features of this boundary 
condition, a good approximation (Greenstein & Hartke 
1983;YP01) is to write it as 

Fr{pb, Qb, B, Tb) « F|| cos^ Qb + F^_ sin^ Ob , (20) 

which expresses Fr, for an arbitrary angle Ob between B 
and the radial direction, in terms of the radial flux for a ra- 
dial field F|| = Fr{pb,OB = 0, B, Tb) and for a meridional 
field F±_ = Fr{pb,OB = 7r/2, 5, Tb). For a dipolar field 65 
is related to the polar angle 9 by tan 9b = 0.5 tan 0. Note 
that when B » 10^^ G one has F± <C i.e., the enve- 
lope is blanketing the star much more strongly around the 
magnetic equator than around the poles. 

2.5. Equation of state and chemical composition 

We will consider a 1.4M0 neutron star whose structure is 
obtained by integrating the Tolman-Oppenheimer-Volkoff 
equation of hydrostatic equilibrium. The core matter is 
described by the equation of state (EOS) calculated by 
Wiringa et al. (1988). We separate the crust from the core 
at the density p = pn = 1.62 x 10^** g cm~'^ (Lorenz et 
al. 1993) and use the EOS of Negele & Vautherin (1973) 
for the inner crust, at p > pdrip = 4.4 x 10"'^'^ g cm~'^, and 
Haensel et al. (1989) for the outer crust. We assume that 
the chemical composition is that of cold catalyzed matter, 
as is the case for these two crustal EOSs. 

The star so produced has a radius Rns = 10-9 km and 
the crust-core boundary is at radius r„ = 0.925i?Ar5 = 
9.33 km. 

2.6. Numerical method 

We solve the heat transport equation in its time- 
dependent form, 

dT 

C^ = -V-F (21) 

starting with constant temperature, T = Tcore- As bound- 
ary condition at the lower boundary the temperature is 
kept fixed to the value Tcore at all times while at the outer 
boundary the radial heat flux is related to the tempera- 
ture at the same point as described above. The code is 
then run until a stationary state is reached. 

For the discretization of the spatial derivatives a stag- 
gered mesh method as described in Stone & Norman 



Geppert, Kiiker & Page: Temperature distribution in magnetized neutron star crusts 



1.0 
0.8 
0.6 



< 

DC 

LU 

I 0.4 

LU 
I- 

0.2 



0.0 



10 



TT 



10^ 10 
DENSITY 



TT 



10 



T4 




10'^ 10 
DENSITY 



10'^ 10' 
DENSITY 



Fig. 5. The crustal temperature T (normalized on Tcore) as a function of the density within the crust at 4 different polar angles: 
6 = 0: dotted-dashed line (almost coincident with the upper limit of the figure), 6 = 30: dashed line, 6 = 60": dotted line, and 
= 90°: full line, for a crustal field of strength Bq = 3 - lO^^G. The three panels correspond, from left to right, to Tcore = 10®K, 
lO'^K, and 10®K, respectively. 



(1992) is used in spherical polar coordinates (r and 9). For 
the integration in time the fully explicit method turned 
out to be prohibitively slow because the diffusion coeffi- 
cients vary strongly in space in the case of strong magnetic 
fields. We therefore use the scheme in an operator-split 
implementation that treats the radial diffusion terms im- 
plicitly while the horizontal diffusion and the dT/dO term 
in Fr remain explicit. This keeps the computational ef- 
fort per time step small as the equations to be solved are 
tridiagonal, but allows for a sufficiently large time step to 
reach a stationary state within several hours of computing 
time. For the evaluation of the transport coefficients and 
the outer boundary condition the temperature resulting 
from the last time step is always used, i.e., only the radial 
derivatives are treated implicitly. 



3. Results 

3.1. Internal temperature 

The importance of the magnetic-field-induced non- 
isothermality of the whole crust is illustrated in Figs. 5,6, 
7, and 8. The temperature profiles of Fig. 5 when com- 
pared with Fig. 6, show clearly the difference between 
a crustal and a core field, the latter inducing temper- 
ature variations in the crust of much less than 1% at 
-Bo = 3 • 10"'^^ G while the former can result in variations 
of a factor two for the same dipolar external field strength 
and Tcore = lO^K, and even much larger at lower Tcore- 

For strong fields, when WbT ^ 1, one has K|| = 
Ko ^ K_L and hence heat flows essentially along the fleld 
lines and, given the large values of kq, no large temper- 
ature gradient can build up along them, as illustrated in 
Fig. 7. Only extremely strong core fields may cause sig- 
nificant deviations from the isothermality of the crust; for 
-Bo = lO^^G we could observe a difference of only 10 % for 
T(pb) between pole and equator. The qualitative difference 
between core and crustal magnetic fields is then easily un- 
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Same as Fig. 5 but for a dipolar core field with Tec 



derstood by observing that field lines are essentially radial 
in most of the crust in the case of a core field while they 
are predominantly meridional for a crustal field inhibit- 
ing radial heat flow in a large part of the crust. We could 
find significant differences to the isothermal crust model 
only if the polar surface field strength exceeds lO^^G. 
Additionally, the core temperature should be smaller than 
lO^K. 

In the case of a star-centered core field. Fig. 6 and 
the right panel of Fig. 7 shows that the stellar equator is 
very slightly warmer than the pole. This is a direct conse- 
quence of the outer boundary condition, i.e., the envelope: 
magnetic-field-induced anisotropies are weak within the 
crust for such fields but are large in the envelope which 
is more insulating around the equator than around the 
poles. 

A 3D representation of the crust temperature is shown 
in Fig. 8: it may be surprising in the sense that heat flows 
into the crust from the core, at a fixed Tcore, and out of the 
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crust at the surface and most of the heat comes out in the 
polar regions which are, first, warmer than the equator and 
where, second, the envelope is less insulating. So heat must 
be flowing from the equatorial regions toward the polar 
regions, i.e., apparently from cold regions toward warmer 
ones ! That this situation cannot violate the second law of 
thermodynamics is built into the heat conductivity tensor 
k which is positive definite and guarantees that F-'VT < 
always. Figure 9 illustrates this situation: at a point in 
the northern hemisphere VT is almost perpendicular to 
B and — (Vr)^ is clearly pointing toward the equator. 
However, since i^y — — K||(VT)||, F± = — k_l(VT)_l, and 
K|| ^ K± we have |i^||| » |i^_L| and the resulting Fg is 
negative, i.e., pointing toward the pole: heat is flowing 
from the equator toward the poles but does it along the 
magnetic field lines from warmer regions toward colder 
ones. Our numerical results show that, e.g., in the northern 
hemisphere, Fg is negative in large regions and is positive 
elsewhere depending on the orientation of B. 




Fig. 7. Representation of both field lines and temperature dis- 
tribution in the crust whose radial scale {r{pn) < r < r{pi,)) 
is stretched by a factor of 5, assuming Bo = 3 ■ 10^^ G and 
Tcore = lO^K. Left panel corresponds to a crustal field, right 
panel to a star-centered core field. Bars show the temperature 
scales in units of Tcoro- 



The results shown in Fig. 5 for the crustal field con- 
figuration with the typical strength Bq = 3 ■ lO^^G con- 
firms the statement that with increasing magnetization 
parameter the anisotropy of the temperature distribution 
within the crust increases too. The magnetization pa- 
rameter increases with increasing magnetic field strength 
and with decreasing crustal temperature, i.e. with de- 
creasing Tcorc, because the relaxation time of electron- 
phonon collisions grows strongly in the course of cooling. 
Therefore, while the temperature profile along the poles 
shows practically no gradient with decreasing Tcorc the 
ratio T{pb,9 = 90°)/T{pb,9 = 0°) decreases from 0.95 to 
0.5 and 0.2 when Tcorc decreases from 10® K to 10^ and 
10^ K, respectively. Also, an increase of the magnetic field 
strength amplifies that difference. Applying the same field 




Fig. 8. 3-D presentation of the temperature distribution in the 
crust for Tcorc = 10^ K and a crustal field with strength Bo = 
3 • 10^'^ G. (corresponding to the right panel of Fig. 5). 




Fig. 9. Illustration of the magnetic-field-induced anisotropy 
(see text for details) 



structure but Bq = lO^^G the temperature ratio becomes 
smaller than 0.1 for Tcorc = lO'^K. 

Note also that the highly unknown parameter of the 
impurity concentration (Q = 0.1 throughout this pa- 
per) affects the relaxation time: the more impurities the 
shorter the relaxation time of electron-impurity collisions. 
Therefore, in a very pure neutron star crust the magnetic 
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field effects onto the crustal temperature distribution will 
be even more pronounced. 

3.2. Surface temperature distribution 

The magnetic field permeating the envelope induces a 
non-uniform surface temperature distribution, mostly due 
to quantizing effects of the field at low densities, even 
in the case of a uniform crustal temperature (Schaaf 
1990a, 1990b; Page 1995). The non-isothermality of the 
crust produced by a crustal magnetic field will result in an 
even more pronounced non-uniformity of the surface tem- 
perature. These effects are shown quantitatively in Fig. 10 
where the two cases of core and crustal fields are com- 
pared. 

The different field structures do not only affect the re- 
lation between polar and equatorial surface temperature 
but also the setup and the extension of the warm polar re- 
gions. In Fig. 10 it is seen that the crustal field can cause 
a much smaller warm polar region than a star-centered 
core field of the same polar surface strength would do. 
While for the core field geometry with Bq = lO^^G and 
Tcore = lO^K the surface temperature at a polar angle of 
31.5° is reduced only to 0.93 of its polar value, for the 
crustal field configuration the corresponding value is 0.82. 
With decreasing T^ore this difference becomes larger: thus 
for Tcoro = lO^K the corresponding values for the core and 
the crustal field are 0.93 and 0.77, respectively. The setup 
of a clearly distinct warm polar region becomes more and 
more pronounced with an increasing magnetization pa- 
rameter (i.e. increasing field strength and/or decreasing 
crustal temperature and/or decreasing impurity concen- 
tration) for a crustal field configuration, while the shape of 
the surface temperature distribution is almost unaffected 
by the magnetization parameter in case of a core magnetic 
field. 

3.3. Luminosity 

Having a non-uniform surface temperature distribution, 
Tg = Ts{9), the effective temperature T^a can be calcu- 
lated, from its definition, as 



B[G] 



(22) 



where R is the circumferential neutron star radius and 
dE = smOdcbAO is the surface area element, so that finally 
the luminosity is given by 



Jo 



L = 2'KasBR^ / T^{e) sin OdO. 



(23) 



The photon luminosities for the neutron star models under 
consideration with different magnetic field structures and 
strengths are listed in Table 1 for a model with Tcore = 
lO^K, Table 2 for Tcoro = lO^K, and Table 3 for Tcore = 
lO^K. 

The almost identical luminosities obtained for the 
isothermal crust and a crust penetrated by a star-centered 



L[erg/s] 
Isothermal 
crust 



L[org/s] 
Star-ccntcrcd 
field 



L[erg/s] 
Crustal 
field 



1.66 • 10^^ 1.64 • 10^^ 1.63 • 10^^ 
10^^-^ 1.85 • 10^^ 1.83 • 10^^ 1.68 • 10^^ 
10" 2.20 -10^^ 2.17-10^^ 1.25-10^^ 

Table 1. Photon luminosities for a neutron star with Tore = 
10*^K. The differences between isothermal crust and star- 
centered field are negligible but significant between them and 
crustal fields larger than lO^^G. 



B[G] 



L[erg/s] 
Isothermal 
crust 



L[erg/s] 
Crustal 
field 



10'^ 

1012.5 
10" 



1.42 • 10 

1.73 ■ 10^ 
5.01 • 10^ 



3D- 



1.31 • 
1.19 ■ 



10 

10^ 



3D- 



1.00 • lO'' 



Table 2. Photon luminosity for a neutron star with Tot 
lO^'K. 



B[G] 



L[erg/s] 
Isothermal 
crust 



L[erg/s] 
Crustal 
field 



T(P 3.36 ■ 10^** 2.82 • lO-" 
4.54 ■ 10^** 2.71 • 10^ 



10^^ 

1& 



0.92 ■ il)-^ 2.5!) • 10^ 



Table 3. Photon luminosity for a neutron star with 
10®K. 



core field (Table 1) reflect the little effect even strong fields 
of that structure have onto the crustal temperature dis- 
tribution. While for a star-centered magnetic field (both 
with isothermal and non-isothermal crust) the photon lu- 
minosity increases with increasing field strength, in neu- 
tron stars possessing a crustal magnetic field above a cer- 
tain strength (« lO^^G) the luminosity is reduced since 
then over the major part of the surface the heat insu- 
lating effect of such a field configuration dominates; its 
^-component causes strong meridional heat fluxes toward 
the polar region whose area, however, becomes smaller 
with increasing fleld strength. This effect impedes the ra- 
dial heat transport strongly, finally less heat can be ir- 
radiated away from the surface and the photon cooling 
process will be decelerated significantly in comparison to 
a non-magnetized neutron star (here Bq < lO^^G) or even 
to a strongly magnetized neutron star with a star-centered 
core field. 

4. Discussion and conclusions 

The insulating effect of a quantizing magnetic field com- 
ponent perpendicular to the radial heat flux is well known 
and its role in the low density layers of the envelope of a 
neutron star has been extensively studied (e.g., Hernquist 
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Fig. 10. The surface temperature Ts as a function of the polar angle 6 and for Tcore = lO^K (left panels), Tcore = lO^K (mid 
panels), or Tcoro = 10*K (right panels) The dashed lines show the surface temperature distribution when an isothermal crust 
is assumed. The full lines represent the surface temperatures when the crust temperature distributions take into account the 
anisotropy of heat transport induced by a crustal magnetic field (the temperature at the crust-core interface being fixed at the 
Tcore)- Almost indistinguishable from the isothermal crust model is the Ts-distribution for a star-centered core field. It is shown 
by dot-dashed lines for Tcore ~ 10*K; for lower Tcore the differences are even smaller. The assumed polar surface field strengths 
-Bo are 10^^ G (upper panels), 3 • 10^^ G (mid panels) and 10^^ G (lower panels). 



1985; Schaaf 1990a; Heyl & Hernquist 1998; YPOl). Here, 
importance of the classical (Larmor) effect for the heat 
transport through the whole crust is demonstrated, pro- 
vided the neutron star possesses a sufficiently strong mag- 
netic field maintained by electric currents circulating in 
the crust but which does not penetrate the core of the 
star. Besides its insulating effect the tangential crustal 
magnetic field creates a meridional heat fiux from the 
equatorial regions towards the high latitude ones. It trans- 
ports the heat, which is dammed in the equatorial region, 
towards the poles where it can be much more easily irra- 
diated away. Eventually this leads to an equatorial belt 
which is much cooler than the poles. These are a direct 
consequence of the confinement of the magnetic field lines 
to the crust which result in the presence of a large merid- 
ional field component in most of the crust. 

Recently, Svidzinsky (2003) argued that the accumu- 
lation of magnetic field lines along the proton supercon- 
ductor at the crust-core boundary, due to the Meissner- 
Ochsenfeld effect, produces an insulating barrier prevent- 
ing heat to fiow between the crust and the core. Our results 
are in the same line of thought but they do not confirm his 
claims of insulation of the crust from the core. The field ge- 
ometries, i.e., the Stoke's stream functions (Eqs. 6, 7), we 
used in our calculations come from field evolution models 
(Page et al 2000) of crustal fields in which the migration, 
and accumulation, of the currents and the field lines to- 
ward the crust-core superconducting boundary was mod- 
eled in detail and, as illustrated in Fig. 7, they allow heat 
diffusion through the crust-core boundary. Nevertheless, 



other field evolution scenarios, e.g. expulsion of the mag- 
netic flux from the core by a proton type I superconductor 
(Link (2003); Buckley et al 2003) may produce a much 
stronger piling up of field lines tangentially to the crust- 
core boundary and result in more efficient thermal insula- 
tion. 

In case the magnetic field is allowed to penetrate the 
core of the star, and assuming a star-centered dipolar ge- 
ometry in the crust, we have shown that the stationary 
thermal state of the crust is very close to isothermality 
thus confirming, for this field geometry, the assumptions of 
the models of surface temperature distribution of magne- 
tized neutron stars (e.g., Grcenstein & Hartke 1983 ; Page 
1995; Page & Sarmiento 1996; Shibanov & Yakovlev 1996) 
which considered only the insulating effect of a quantizing 
field in the envelope and assumed the rest of the crust to 
be isothermal. 

The obvious difference between the isothermal crust 
and our crustal field results is demonstrated in Fig. 5 and 
6. For a crustal field these differences cannot be neglected 
for highly magnetized neutron stars. 

Since the warm polar region in case of a strong crustal 
field is much smaller than for star-centered and/or weak 
magnetic fields (see Fig. 10) this may open a new way to 
distinguish between crustal and core magnetic fields: A 
strong crustal magnetic field implies a smaller effective 
area for thermally emitting cooling neutron stars. This 
has consequences which can potentially be observed in 
X-ray in cases where fits of the thermal component of 
the spectrum of a cooling neutron star result in effective 
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emitting radii i?eff which are significantly smaller than 
the expected radius of a neutron star. For example, the 
"Three Musqueteers" (PSR 0656+14, PSR 1055-52 and 
Geminga: Triimper & Becker 1998) as well as RX J185635- 
3754 (Pons et al 2002) all have iieff ~ 5 - 7 km when their 
thermal spectra arc fitted with blackbody spectra. If these 
radii would coincide with the radius of the neutron star, 
the equation of state describing the state of the core mat- 
ter would have to be extremely soft. A relatively small 
warm polar region, created by a strong crustal field and 
emitting almost all the thermal radiation would be a rea- 
sonable explanation for such small i?off- 

The differences in the photon luminosities for a star- 
centered or a crustal field will also affect the long term 
cooling of neutron stars. Future cooling calculations as 
well as the comparison of X-ray spectra and pulse profiles 
with model calculations assuming different field structures 
will open a way to discriminate the two basic scenarios: 
crustal or core magnetic field. Finally we mention the con- 
sequences the non-isothermality of the crust may have for 
the crustal field itself. Since the electric conductivity in 
the hot polar region is much smaller than in the equa- 
torial layer, the field decay will be affected too and may 
cause differences in the field structure. This "back reac- 
tion" of the field onto its own decay via a field driven 
non-spherical symmetric crustal temperature distribution 
is subject of future investigations. 
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